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ABSTRACT 

We have carried-out multi-configuration Breit-Pauli AUTOSTRUCTURE calculations for the dielec- 
tronic recombination (DR) of Fe*+ - Fe^^+ ions. We obtain total DR rate coefficients for the initial 
ground-level which are an order of magnitude larger than those corresponding to radiative recombi- 
nation (RR), at temperatures where Fe 3p'' (g = 2 — 6) ions are abundant in photoionized plasmas. 
The resultant total (DR+RR) rate coefficients are then an order of magnitude larger than those cur- 
rently in use by photoionized plasma modeling codes such as CLOUDY, ion and xstar. These rate 
coefficients, together with our previous results for q — Q and 1, are critical for determining the ion- 
ization balance of the M-shell Fe ions which give rise to the prominent unresolved-transition-array 
X-ray absorption feature found in the spectrum of many active galactic nuclei. This feature is poorly 
described by cloudy and ION, necessitating an ad hoc modification to the low-temperature DR rate 
coefficients. Such modifications are no longer necessary and a rigorous approach to such modeling can 
now take place using these data. 

Subject headings: atomic data - atomic processes - galaxies: active - galaxies: nuclei - X-rays: galaxies 



L INTRODUCTION 

The iron M-shell ions (predominantly Fe^"'' - Fe"'^'^+) 
give rise to strong X-ray absorption lines due to n = 
2 — 3 electronic inner-shell transitions. These are seen 
in the spectra of active galactic nu clei (AGNs) obser ved 
by Chandra and XMM-Newton l)Sako et alJ l2001|) as 
an unresolved-transition-array (UTA). The shape of the 
UTA feature can b e used as a powerful diagnostic — see 
iBehar et al.l l)200H) for a detailed discussion, and whose 
atomic d atabase for desc ribing it ha s been upda t ed re- 
cently bv lGu et all (2006). However, Nctzcr et al. ( 2003') 
have pointed- out problems in modeling this shape an d so 
iNetzed lj2004D . using ION, and lKraemer et all lj2004|l . us- 
ing CLOUDY, have suggested increasing the magnitude of 
the recombination rate coefficients for Fe ions, particu- 
larly the 3p'' [q = 1—6) ions Fe^'^+ - Fe^"'', by postulating 
ad hoc low-temperature dielectronic recombination (DR) 
rate coefficients. The net effect is to shift the ionization 
balance towards the neutral end. This brings the mod- 
eling results towards accord with observation, but it is 
not rigorous since there is little to constrain the DR rate 
coefficients which they use. 

The extant recombination data is effectively all 
radiative (RR) at photoionized plasma tempera- 
tures ^. The recommende d DR data for these 
Fe ions ijArnaud fc Ravmondl '1992) is based upon 
th e high-temper a ture e lectron-coUisional plasma results 
of I.Tacobs et al.l l|1977 ^. which give essentially zero- 
contribution from DR at photoionized plasma temper- 
atures. There is much evidence from the Fe L-shell, 
both experimental and theoretical (see, e.g.. lSavin et al.l 
I2OO 6'). that a significant contribution can be expected 
from DR at low temperatures, due principally to 'non- 
dipole' core-excitations which were not (needed to be) 

^ When we refer to photoionized or electron-collisional plasma 
temperatures we mean the tem peratures at which Fe 3p'^ ion s are 
abundant in such plasmas — seelKallman and Bautistal 12001]) (xs- 
tar) and lMazzotta et alj 119981) . respectively. 



considered bv lJacobs et a l.l (1197^ 

Recently, measurements by Sc hmidt et all II2006D ^ on 

Fe^'^^ (g = 1) at the Heidelberg heavy-ion test storage 
ring have found the contribution from DR to the total 
recombination rate coefficient to be an order of magni- 
tude larger than that due to RR, at photoionized plasma 
temperatures. We have carried-out a detail ed theo- 
ret ical analysis of these experimental results ijBadnell 
l2006ai) . We found some disagreements, especially at en- 
ergies which affect the DR rate coefficient at photoion- 
ized plasma temperatures. Nevertheless, our DR rate 
coefficient for Fe"'^^+ is also an order of magnitude larger 
than our RR one, but it is up to a third smaller than 
the experimentally-based one over lO'' — 10^ K. We now 
report-on the results of calculations for the computation- 
ally demanding Fe 3p' (g = 2 — 6) ions. 

The remainder of the paper is organized as follows: 
in Section 2 we describe our methodology, in Section 3 
we present our results and make comparisons with the 
results of earlier works, and then make some concluding 
remarks. 

2. METHODOLOGY 

Our recent comprehensive study l|Badnell2006a|) of the 
DR of Fe^'^+ (g — 1), including detailed comparisons 
of our theoretical cross sections with th ose of th e high- 
energy resolution measurements of Schmi dt et al.l (j2006) , 
provides us with a guide for our approach to the 3p'' 
(g = 2 — 6) ions. These are both computationally more 
demanding and there are no (published) experimental 
results for them. We give the essentials of the new work 
below and refer the reader to Badnell ( 20 06aj) for a more 
detailed exposition of the atomic physics. But, we note 
in particular that a level-resolved treatment is critical 
for the determination of accurate DR rate coefficients 
at photoionized plasma temperatures, but the results of 

2 ISchmidt et"an l2()(13) also give an up-to-date list of references 
of observations of the AGN X-ray absorption feature. 
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such calculations are sparse for the M-shell. 

2.1. Theory 

The total dielectronic recombination rate coefficient, 
aJ^^(T), from an initial state v of an A^-elec tron ion is 
given, at a temperature T, by l)Burgesglll964j) 



,DR 



(r) = 



3/2 



/47ra|iH 



^ 2w 



(1) 



where o^j is the statistical weight of the [N + l)-electron 
doubly-excited resonance state j, oJi, is the statistical 
weight of the initial state and the autoionization [A'^) 
and radiative {A^) rates are in inverse seconds. Here, 
is the energy of the continuum electron (with orbital an- 
gular momentum /), which is fixed by the position of the 
resonances, and /h is the ionization potential energy of 
the hydrogen atom (both in the same units of energy), 
A;b is the Boltzman constant, T the electron temperature 

and (47rag)3/2 = 6.6011 x IQ-^^ cm^. 

We have used autostructure (jBadnelM |19.86D to 
carry-out multi-configuration Breit-Pauli calculations of 
all of the necessary rates and energies, as detailed next. 

2.2. The Fe 3p« target 

We describe the A^-electron target by the following con- 
figurations (assuming a closed Ne-like core): 

1: Ss^Sp?, 2: SsSp^+i , 3 : Ss^Bp^-iSd, 

4 : 3p9+2 , 5 : SsSp'Sd , 6 : Ss^Sp^-^sdS 

7: Sp-J+^Sd. 

(If a given value of q results in an occupation number 
< or > 6 then that 'configuration' does not exist for 
said ion.) This target expansion allows for both 3s and 
3p An — sub-shell promotions from the ground config- 
uration, as well as including important interacting con- 
figurations. We denote it as '7CF'. We also investigated 
the accuracy of this representation for DR by carrying- 
out a further calculation for each ion which included the 
additional configuration interaction due to 

8 : 3s3p«-i3d2 . 

We denote the combined set as '8CF'. 

The contribution from higher energy (An = 1) promo- 
tions has no effect at photoionized plasma temperatures, 
while in electron-collisional plasmas we expect such con- 
tributions to be less than 10% of the An = 0, follow- 
ing our results for q = 1 (BadncU 2006a). (The relative 
importance of An — 1 contributions decreases as q in- 
creases. This is because the inner-shell 2 — 3 contribution 
is further suppressed by additional core re-arrangement 
autoionizing transitions and the 3 — 4 outer-shell contri- 
bution peaks at a temperature which is closer to that of 
the An = peak.) 

All relevant {N + l)-electron autoionization and ra- 
diative rates are then determined, for the given target 



expansion, for all Rydberg states up to n = 1000 and 
/ = 9 — 13, depending on q, and the total DR rate coef- 
ficient is then determined according to equation ( 1). 

In addition, we use observed target energies 
lv3D wherever possible. This minimizes the sensitivity to 
the Maxwellian exponential factor at low temperatures, 
which is critical for application to photoionized plasmas. 
Since the observed energies are incomplete, we adopt the 
expedient strategy of using calculated level splittings to 
adjust the position of missing levels of a term, relative 
to an observed one. Similarly, we adjust (unobserved) 
terms to maintain the splitting with observed ones of 
the same symmetry. Finally, for higher energy config- 
urations (typically, configuration number 4 and above) 
where there are no observed level energies, we adjust the 
entire configuration position by consideration of the av- 
erage shift of lower- lying configurations. 

Radiative recombination may be expected to be im- 
portant at photoionized plasma temperatures and so we 
have also calculated total RR rate coefficients, aJ^^(T), 
with AUTOSTRUCTURE, following Badnell ( 2006bfl . 

2.3. Fits 

For ease of use, we fit our total recombination rate 
coefficients to the usual functional forms: for DR, 



,DR 



(T) = T-3/2^c,exp(-ii;,/T) 



(2) 



where the Ei are in the units of temperature, T (K), and 
the units of q are then cm^s^^ K^^^; for RR, 



a^'\T) = A 



(3) 



where, for greater accuracy, B may be replaced as 

B ^ B + Cexp{~T2/T). (4) 

Here, 7o_i^2 a-i'e in units of temperature (K), the units of 
A are cm^ s~^, while B and C are dimensionless. 

3. RESULTS 

In Figures 1 - 5 we present our recombination rate co- 
efficients for the initial ground-level of Fe"'^^+ - Fe^+ {q — 
2 — 6) and compare them wi th the D R rate coefficients 
re commended b y Arnaud & Ravmondl f 1992(1 . The data 
of lArnaud fc Rav mond (. 1993) are based principally upon 
the results of iJacobs et all l|1977|) . but include an esti- 
mate of the contribution from 2p — 3d inner-shell transi- 
tions as well — thes e were not in c lude d by Jacobs ct aQ 
l(1977|) . The work of lJacobs et alJ lfl977f) was for applica- 
tion to high-temperature electron-collisional plasmas and 
their results have little in the way of a low-temperature 
contribution — a modest one may just be seen for g = 4 
and 5 at around 10^ K. 

The main result, common to all ions, is the order of 
magnitude difference between the DR and RR rate co- 
efficients at photoionized plasma temper atures. This is 
similar to the situation found for q = 1 iSchmidt et alj 
2006: Badnell 2006a). Thus, the ad hoc modifications 
proposed bv lNetzeii (j20041|l and bv lKraemer et all l|200l) 
were conservative — increasing the total recombination 
rate coefficient by factors of only 2-4. These new rate 
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coefficients (plus q = 0,1 data) can be expected to have 
a significant effect upon the ionization balance of Fe^+ - 
Fe^'^"'" in photoionized plasmas and, in turn, the model- 
ing of the UTA X-ray absorption feature in AGNs. On 
comparing our 7CF and 8CF results, we note also that 
sensitivity to near threshold resonance positions at low 
temperatures does not appear to become an issue until 
< 10'* K, i.e., below the main temperatures of interest 
for photoionized plasmas. 

At electron-coUisional plasma temperatures we find 
reasonable agreement with the recommended data of 
lArnaud fc Ravmondl l|1992f) . The strength of the contri- 
bution from 'low-temperature' resonances extends to en- 
hancing the high-temperature peak in many cases. Even 
at 10^ K, which is far-off equilibrium, our present re- 
sults are still (a little) large r than those recommended by 
lArnaud fc Ravmondl l)1992() . except for the case of g = 2. 
Even in this case {q = 2) we expect the An = 1 contri- 
bution to result-in no more than a 10% increase at a few 
times 10^ K, rising to at most 20% by 10'' K, and less 
for lower charg e ions, based-upon our results for q = I 
(|Badnel]|l2006a|) and our comments in Section 2.2. Thus 
the present results can, and should, be used for modeling 
electron-coUisional plasmas as well. 

In Tables 1 and 2 we present our fit coefhcients for these 
DR and RR rate coefficients, respectively, as defined by 
Equations (2) an d (3). For conv enience, we include the 
results for (J = 1 l)Badnellll200"6a|) and g = l|Altun et all 



12006ft as well. Results f or all hig her charge states may 
be found online ( Badnclil l2fl06cldfl' following the work of 
iBadnell et all l|2005) and iBadnelll l|2006bO for DR and 
RR, respectively. 



4. CONCLUDING REMARKS 

We have reported new DR rate coefficients for Fe 3p* 
(q = 2—6) ions which give rise to total recombination rate 
coefficients which are an order of magnitude larger, at 
photoionized p lasma temperatures, than t hose currently 
recommended ()Arnaiid fc R.avmonc]llT992h and routinely 
used by modeling codes such as CLOUDY, ion and xstar. 
These new rate coefficients can be expected to change 
significantly the ionization balance of the Fe M-shell ions 
which give rise to the important UTA X-ray absorption 
feature seen in the spectra of AGNs observed by Chandra 
and XMM-Newton, and to enable rigorous modeling of 
it to be carried-out now using these data. 

Similarly, large low-temperature DR contributions can 
be expected for 3p'' ions of other elements of astrophys- 
ical importance, e.g., Si, S, Ar, Ni. A move into the 3d 
sub-shell of lower-charge Fe ions, using a level-resolved 
approach, is also desirable. 
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TABLE 1 

DR FITTING COEFFICIENTS, Cf (cm^S~-'K^/^) AND Ef(K), FOR THE INITIAL GROUND-LEVEL OF Fe 3p« ((J = — 6) IONS. 



q 


Cl 


C2 


C3 


C4 


C5 


C6 


C7 


C8 





5.636(-4)=' 


7.390(-3) 


3.635(-2) 


1.693(-1) 


3.315(-2) 


2.288(-l) 


7.316(-2) 




1 


1.090(-3) 


7.801(-3) 


1.132(-2) 


4.740(-2) 


1.990(-1) 


3.379(-2) 


1.140(-1) 


1.250(-1) 


2 


3.266(-3) 


7.637(-3) 


1.005(-2) 


2.527(-2) 


6.389(-2) 


1.564(-1) 






3 


1.074(-3) 


6.080(-3) 


1.887(-2) 


2.540(-2) 


7.580(-2) 


2.773(-l) 






4 


9.073(-4) 


3.777(-3) 


1.027(-2) 


3.321(-2) 


8.529(-2) 


2.778(-l) 






5 


5.335(-4) 


1.827(-3) 


4.851(-3) 


2.710(-2) 


8.226(-2) 


3.147(-1) 






6 


7.421(-4) 


2.526(-3) 


4.605(-3) 


1.489(-2) 


5.891(-2) 


2.318(-1) 






? 


El 


E2 


Es 


E4 


E5 


Ee 


Er 


Es 





3.628(3) 


2.432(4) 


1.226(5) 


4.351(5) 


1.411(6) 


6.589(6) 


1.030(7) 




1 


1.246(3) 


1.063(4) 


4.719(4) 


1.952(5) 


5.637(5) 


2.248(6) 


7.202(6) 


3.999(9) 


2 


1.242(3) 


1.001(4) 


4.466(4) 


1.497(5) 


3.919(5) 


6.853(5) 






3 


1.387(3) 


1.048(4) 


3.955(4) 


1.461(5) 


4.010(5) 


7.208(5) 






4 


1.525(3) 


1.071(4) 


4.033(4) 


1.564(5) 


4.196(5) 


7.580(5) 






5 


2.032(3) 


1.018(4) 


4.638(4) 


1.698(5) 


4.499(5) 


7.880(5) 






6 


3.468(3) 


1.353(4) 


3.690(4) 


1.957(5) 


4.630(5) 


8.202(5) 






'^Note, (n) denotes xlO". 






















TABLE 2 














RR FITTING COEFFICIENTS FOR THE INITIAL GROUND-LEVEL OF Fe 3p' (q = — 


6) IONS. 




? 


A (cm3 




B 


To(K) 


Ti(K) 




C 


T2(K) 





1.179(- 


^9) 


0.7096 


4.508(2) 


3.393(7) 




0.0154 


3.977(6) 


1 


1.050(- 


-9) 


0.6939 


4.568(2) 


3.987(7) 




0.0066 


5.451(5) 


2 


9.832(- 


-10) 


0.7146 


3.597(2) 


3.808(7) 




0.0045 


3.952(5) 


3 


8.303(- 


-10) 


0.7156 


3.531(2) 


3.554(7) 




0.0132 


2.951(5) 


4 


1.052(- 


-9) 


0.7370 


1.639(2) 


2.924(7) 




0.0224 


4.291(5) 


5 


1.338(- 


-9) 


0.7495 


7.242(1) 


2.453(7) 




0.0404 


4.199(5) 


6 


1.263(- 


-9) 


0.7532 


5.209(1) 


2.169(7) 




0.0421 


2.917(5) 
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Fig. 1. — Total ground-level rate coefBcients for Fe^-^+ (q = 2). Solid curve, DR (7CF); dotted curve, PR (8C F); long-dashed curve, RR; 
all present AUTOSTRUCTURE results. Short-dashed curve, recommended DR data of lArnaud Hz Ravmond' fl993). PP and CP denote typical 
photo ionized and electron-collisional plasma temperature ranges, respectively, for Fe^'^"'" IRa llman an d Bautistal2001l) and JMazzotta et al] 
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Fig. 2.— As Figure 1, but for Fe"+ (q = 3). 
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Temperature (K) 

Fig. 3.— As Figure 1, but for Fe^°+ {q = 4). 
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Fig. 4.— As Figure 1, but for Fe^+ (g = 5). 
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Fig. 5.— As Figure 1, but for Fe^+ {q = 6). 



